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Abstract: We study thermodynamics of free SU{N) gauge theory with a large 
number of colours and flavours on a three-sphere, in the presence of a baryon number 
chemical potential. Reducing the system to a holomorphic large- matrix integral, 
paying specific attention to theories with scalar flavours (squarks), we identify novel 
third-order deconfining phase transitions as a function of the chemical potential. 
These transitions in the complex large- saddle point configurations are interpreted 
as "melting" of baryons into (s) quarks. They are triggered by the exponentially 
large (~ e^) degeneracy of light baryon-like states, which include ordinary baryons, 
adjoint-baryons and baryons made from different spherical harmonics of flavour fields 
on the three- sphere. The phase diagram of theories with scalar flavours terminates 
at a phase boundary where baryon number diverges, representing the onset of Bose 
condensation of squarks. 
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1. Introduction and summary 

Unravelling the thermodynamic phase structure of gauge theories such as quantum 
chromodynamics (QCD) poses several challenging and interesting questions. The 
behaviour of cold, dense matter in QCD is particularly elusive due to the so-called 
sign problem, and because low temperature, flnite density transitions are expected to 
occur in the regime of strong gauge coupling [Q. The large- expansion promises to 
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be a useful means of understanding certain aspects of the relevant physics. In recent 
years, the 't Hooft and Veneziano large- iV limits p|, ^ have been applied to argue 
the existence of a low temperature quarkyonic phase Q where baryons condense. 
A second line of investigation of large- limits of QCD-like theories, involves the use 
of tools provided by gauge/gravity dualities to describe strongly interacting 
theories with a large number of flavours ||^, ^. 

Drawing motivation from these research directions, in this paper we study the 
phase diagram of SU{N) gauge theories with both fundamental scalars and fermions 
in the Veneziano large- limit, on a spatial three-sphere, and with the gauge coupling 
set to zero. For asymptotically free theories, this is a natural limit to consider when 
the radius of the sphere is much smaller than the dynamical length scale at which 
gauge interactions become strong. Despite being in a compact space, a non-trivial 
phase structure is possible because of the large- A^ limit which plays the role of the 
thermodynamic limit. 

It is now well appreciated, following the works of pU] , |n[], that even free large- A^ 
theories on compact spaces can have a non-trivial phase structure which may, in 
certain cases, be continued through to the interacting situation. For theories with 
adjoint matter, most famously for the free JV = 4 supersymmetric Yang-Mills (SYM) 
theory on S^, this study reveals a Hagedorn density of states leading to a first order 
deconfinement phase transition with increasing temperature^. The deconfinement 
transition survives the continuation to strong 't Hooft coupling wherein the gravity 
dual displays a Hawking-Page transition from thermal AdS space to an AdS black 
hole geometry [|T^. The qualitative agreement remarkably persists upon inclusion of 
chemical potentials for R-symmetry charges in A/" = 4 SYM [1^, |T6|. 

We consider models containing both fundamental scalars ( "squarks" ) and fermions 
( "quarks" ) as this is the generic matter content of supersymmetric theories that arise 



in holographic descriptions of gauge theories with fundamental flavour fields (e.g. [[Tf 
and the D3-D7 system 0, [l^). Thermal transitions (at zero density) in such setups 
at zero and weak couplings were considered in ^ ^ The crucial new fea- 
ture upon introducing a chemical potential for (s) quark number or baryon number, is 
that the effective action becomes complex due to the breaking of charge conjugation 
symmetry. The thermodynamics of gauge theory on is encoded in the effective 
action for the Polyakov loop matrix U, or the holonomy of the gauge field around the 
Euclidean thermal circle. This effective theory on a compact space yields a matrix 
model for U. The effective action being complex, it turns out that the matrix integral 
must be interpreted in a certain holomorphic sense so that the integral picks up, and 
is dominated by, complex saddle point configurations at large- A^ (see e.g. p3|)- This 



^ While it is not known whether this first order transition in Af = i SYM continues to exist for 
small non-zero 't Hooft coupling, the corresponding weak coupling calculation was carried out for 
pure SU (N) Yang- Mills theory (at large N) and the first order deconfinement transition found to 
persist p^. 
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was first shown in where theories with fundamental fermions were studied (see 
also 



^ In this paper we will use the same techniques revealing a wider range 
of phenomena that occur when fundamental scalars are included. We further clarify 



the physical reasons behind these and also some of the results found in p4|. 



Bose condensation 




Figure 1: General features of the phase diagram. Thick black lines correspond to third 
order phase transitions. Sawtooth features are associated to the filling of low lying Fermi 
levels. The first such feature appears on either side of N = niBp, the mass of the lightest 
baryon made from fermions (quarks). The value N ^ = niB, the mass of the lightest baryon 
constituted of scalar modes, represents the onset of Bose condensation of squarks and is 
the boundary of the phase diagram for the free theory. With interaction potentials typical 
of supersymmetric models, the theory for Nfi > niB will enter an unstable Higgs phase 
(shaded region), along with a locally stable deconfined phase (thatched region). 



Our results are summarized in Figure(|l]), which provides a general picture of 
the phase diagram for these theories. The theory at /x = 0, with increasing tempera- 
ture, experiences a third order Gross- Witten-Wadia (GWW) ^ phase transition 
characterized by the formation of a gap in the large- eigenvalue distribution of the 



unitary matrix U. The non-analyticity at /i = was first noted in [|T9|. Although 
there is no natural order parameter for confinement at finite in the presence of 
fundamental matter fields, non-analytic behaviour in the Polyakov loops occurs in 
the large- A^ Veneziano limit and can be interpreted as deconfinement of mesons to 
quarks. 

The new features in the phase diagram, at low temperature and non-zero chemi- 
cal potential, can be explained as follows. By Gauss's law on S^, the grand canonical 
partition function must be thought of as a sum over gauge-singlet states only. A pos- 
itive chemical potential fi > essentially favours baryons over anti-baryons. There 
are primarily three distinct types of baryonic states whose masses and degeneracies 
determine the nature of the phase diagram: 



• Ordinary baryons Bp constituted of fundamental fermions or quarks. 
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Baryon operators Bs, composed of fundamental scalars and in particular, var- 
ious spherical harmonics of squark fields on S^. 



• The third type of states which will be relevant for our discussion are the 
"adjoint-baryons", i?Ad- Adjoint baryon operators involve any number of in- 
sertions of adjoint fields into an ordinary baryon operator. The adjoint fields 
could be harmonics of the gauge field on or the modes of any other adjoint 
matter field in the theory. 



At large N, the degeneracy of all such states grows exponentially with A^, as we 
will see below (for related discussions see This exponential growth of baryon 

states has a very important effect: although baryons are exponentially suppressed 
due to their masses being 0{N), their exponentially large degeneracies result in 
a competition between energy and entropy, wherein the latter can overwhelm the 
former in the grand canonical ensemble. Consequently, for any fixed low temperature 
we may argue the existence of a phase boundary at a value of /x determined by the 
mass and the degeneracy of the lightest baryon state on S^. 

The expectations above are confirmed and complemented beautifully by the ex- 
plicit behaviour of the complex saddle points of the large- A^ matrix model for U. As 
H approaches the mass of the lightest baryon state, and precisely when the entropy 
factor overwhelms the Boltzmann suppression, the distribution of eigenvalues of U 
along a closed contour in the complex plane, experiences a "gapping" transition of 
the GWW type. This is accompanied by a smooth increase in baryon number, and 
the phase transition is found to be of third order. The scaling of the grand potential 
with the number of flavours indicates that the confined phase should be attributed 
to an ensemble of gauge-singlet mesons and in this sense the quarks are confined. 
Across the gapping transitions we find that gauge-invariant order parameters holo- 
morphic in U, namely (;^Tr [/"), acquire non-zero values. Non- vanishing expectation 



values for Polyakov loop observables are expected in a deconfined phase [^]. For 
this reason we term the non-analyticities as deconfinement transitions across which 
baryons "melt". The anti-holomorphic observables, (;^Tr W^) remain non- vanishing 
in all phases reflecting the asymmetry under charge conjugation in the presence of a 
(positive) chemical potential. 

If the lightest state(s) happen to be ordinary baryons Bp, composed of light 
quark modes, two successive third order transitions occur on either side of the baryon 
mass at Nfi = T T In T){Bp). Here D is the degeneracy of the light baryon 
state. The first is a gapping transition followed immediately by a transition back to 
the ungapped or confined phase, with increasing fi. Across the two transitions, the 
lowest lying Fermi levels get filled and the baryon number increases by a discrete 
step. In between two such transitions, a 'spike' of deconfined phase penetrates all 
the way down to zero temperature. These are the features shown in Figure(0). 
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In theories with squarks, an additional feature appears in the phase diagram. 
For baryons made from scalar flavours, complete antisymmetry of colour indices 
means that the lightest baryon cannot always be constituted solely of the lightest 
squark harmonics on S*^ (the zero modes). When the number of squark flavours is 
smaller than the number of colours, non-trivial baryon operators Bs can exist only 
if higher harmonics of the squark fields are also included. Alternatively, one may 
form non- vanishing baryon operators -BacI by insertions of adjoint fields. In all cases, 
our matrix model calculation explicitly yields a third order deconfining transition 
at Nfi = rriB — Tin 2)(-B), where B is the lightest baryon state composed of scalar 
flavours and D{B), the associated degeneracy. Beyond this transition, we find that 
the baryon number grows without bound and eventually diverges at Nfi = uib ■ We 
identify this as the Bose condensation of squarks/baryons and the phase diagram 
cannot be continued past this boundary in the free theory. 

The onset of Bose condensation of squarks is accompanied by the collapse of the 
distribution of eigenvalues of f/ to a point, so that U = 1. This marks the entry 
of the field theory into a "Higgs phase" where gauge-invariant operators carrying 
(positive) baryon number could acquire vacuum expectation values (VEVs). It is 
not possible to see this in the non-interacting theory. If interactions, such as those 
in supersymmetric QCD-like theories^, were to be switched on, flat directions in 
the scalar potential would lead to runaway, unstable VEVs for large enough baryon 
number chemical potential 0, 0, 0. At any finite temperature however, there 
would be a metastable phase (thatched region in Figure (|1|)) of deconfined (s) quarks 
at the origin of the Higgs branch |T^, 0. In the fi — T plane, the width of the 
metastable region, at weak coupling, is determined by the thermal mass of the squark 
modes. 

Our motivation for exploring the large- phase structure above is to glean infor- 
mation on the possible thermodynamic behaviour of such theories, which could be 
used as a point of reference for investigating similar theories at strong coupling within 
the framework of gauge/string dualities While it may be unreasonable to ex- 



pect a putative strongly coupled gravity dual to reproduce the features in Figure(|l]), 
certain aspects are possibly robust enough to exist at any coupling. The exponential 
(in A^) growth of baryonic states appears unlikely to be changed by interactions; 
baryon-baryon and meson-baryon interactions are intrinsically large, involving pow- 
ers of A^ but not exponentials The cold deconfined phase at finite density in 
the phase diagram suggests a dual, extremal black hole geometry at strong coupling. 

The organization of the paper is as follows: In Section 2, we introduce notation 
and summarize the matter content of the theories in question. Section 3 is devoted 
to the basic features/transitions in the holomorphic large- A^ matrix model, obtained 



^Note that the formulation of rigidly supersymmetric, non-conformal gauge theories on curved 
spaces, such as S*^ x M is subtle as it requires the introduction of new terms in the action 04, Ra, am 
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from the theory truncated to hghtest modes. We also show how the observed non- 
analyticities coincide with the expectations from counting of baryon operators. In 
Section 4, we show that in theories with a small number of squark flavours, it becomes 
necessary to include higher harmonics and heavier modes to obtain the correct phase 
diagram. In Section 5 we review the fi = 0, finite T behaviour of the partition 
function. We end with Section 6, a summary of the phase diagram and possible 
lessons for related systems at strong coupling. 



2. Free gauge theory on 

To evaluate the grand canonical partition function of free gauge theories, we use 
the approach of |rT|. We consider SU{N) gauge theory with upN Dirac fermions 
("quarks"), and ugN complex scalars ("squarks"), all transforming in the fundamen- 
tal representation of the gauge group. The number of flavours is chosen to naturally 



scale with A^. The theory with z/5 = has already been studied in ||2J]. Below we 
extend this to include fundamental scalars, and further provide a natural interpre- 
tation of the ensuing phase diagram as a function of temperature and baryonic or 
quark chemical potential. The thermal partition function of the gauge theory on 5*^ 
is computed by formulating the field theory in Euclidean signature, on x with 
the 5*^ circumference 13 = ^, playing the role of the inverse temperature. 

For the free field theory on S^, the partition sum follows from one loop fluctua- 
tion determinants for the Kaluza-Klein harmonics of all fields on subject to the 
constraint from Gauss' law. This latter constraint is necessary since the sum must 
be over physical, gauge-invariant states, and in practice, it is elegantly implemented 
by introducing a background thermal Wilson line for the zero mode of the gauge field 
on S^: 



U = 6'" 

a = 



Vol(S3 



/ Ao = diag(^i,^2,...,^7v), V^a = 0mod27r.(2.1) 

SU{N) gauge transformations can be used to put this Wilson line - also known as 
the Polyakov loop, in a diagonal form as above. The problem of calculating the 
partition function is then reduced to determining the quantum effective action for 
the "light" zero modes above, after integrating out all Kaluza-Klein harmonics on 
as quadratic fluctuations. 

We can consider the general situation where our 3+1 dimensional gauge theory 
has both adjoint and fundamental matter fields ^. Let ns and np be the number of 

^This is a subtle issue once interactions are switched on and the gauge couphng acquires scale 
dependence. For an asymptotically free theory we must ensure that the radius of the 5''^ is sufficiently 
small compared to the relevant strong coupling scale. If the theory is infrared (IR) free, then the 

must be taken to be large, and the flavour masses small compared to the scale of any putative 
Landau pole. 
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adjoint (real) scalars and chiral fermions respectively, 



0,, /=l,2,...n5, (2.2) 
ipA, A = 1,2,... Up. 

For simplicity we take the adjoint fermions and scalars to have vanishing bare mass, 
with the scalars conformally coupled to the curvature of the three-sphere. 

In addition, we take 2A^5 complex scalars ("squarks") in the fundamental rep- 
resentation, each with mass m^, and Np fundamental Dirac fermions ("quarks") of 
mass uiF, 

[q, , q,) : {N, N) of SU{N) i = l,2,...Ns (2.3) 
(V^, , ^i) : (N, N) of SUiN) t = l,2,...NF. 

Our labelling of the flavour fields is motivated by supersymmetric theories with fun- 
damental matter (SUSY QCD), although the analysis we present is general without 
any reference to supersymmetric theories. We also take all scalar flavours ( "squarks" ) 
to be non-minimally coupled to the curvature of 5*^, which ensures that the fluctua- 
tions on 5''^ have a gap even when the flavour masses are set to zero. In a conformal 
field theory (as is the case for certain supersymmetric examples), the scalars would 
be massless but conformally coupled to the curvature, ensuring a gap. 

We define a U{1)b global symmetry, under which the fundamental fermions and 
scalars {qi,'ipi) carry a charge +1, whilst the anti-fundamental fields {qi,4'i) have 
charge —1. We will refer to this as "quark number" or "squark number". The 
number of baryons is related to it by an extra factor of Introducing a chemical 
potential for the squark or quark number, the effective action for the theory on a 



spatial takes the form of a unitary matrix model (see e.g. JIT], |T4|, |T5|, ^ ^), 

^53x51 = J [dU] exp{-S,s[U]) . (2.4) 

The effective action Sqs consists of two separate contributions, one from the adjoint 
sector and the other from the flavour sector, taking the general form 

•S'cff = •S'adj + "Sfund (2.5) 

oo 

= -J^i-^^ Tr t/" Tr f/t" + j; (Tr t/" e''^^ + Tr f/^" e""^'^)] . 

n=l 

The adjoint fields are uncharged under U{1)b and the ensuing contributions are 
proportional to TtW^TtU^^. These terms are invariant under large gauge transfor- 
mations that "wind" around the thermal up to an element of Z^y. Fundamental 
flavour fields break the Z^r invariance explicitly and give rise to single-trace terms 
of the type Trf/"' and Trf/"''". The most important point to note is that, due to 
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the non-zero baryon number chemical potential, the effective action is complex. The 
coefficients An and can be explicitly written in terms of single particle partition 
functions Zp,zs and zy., for the fermion, scalar and vector harmonics respectively, 
on the three- sphere. 

An = -[zv{x^) + nszs{x^) - {-irnpZFix'')] (2.6) 
n 

X = zv{x) = — zs{x) = zf{x) = 



:i-x)3 ' (l-a;)3 ' (l-x)3 ■ 

Here R is the radius of 5''^. The degeneracies and energies eg of the different modes, 
which lead to each of the partition sums above are summarized in Table (|l]) in the 
Appendix. For the flavour modes the coefficients J-'n are given as 

J'n = I [Ns Zs (n|, msR) - Np (-1)" Zp (n|, tufR) ] . (2.7) 

As before Zs and Zp are partition functions of scalar and fermion harmonics on 
with non- vanishing bare masses ms and rriF respectively: 



(2. 



Zs{ji,msR) = J2 + 1)' (-f + 2) + S + m|i?2^ 
Z^(|, m^/?) = f^2£(£+ l)exp + |)' + m|i?2^ ^ 

where S denotes the coupling to the Ricci scalar of S^. For conformal coupling H = 1, 
while minimal coupling corresponds to S = 0. 

The effective action in the adjoint sector scales as N'^, for large N (taking Ug ~ 
Up ~ C^(l)) as each term in the action is a double trace operator, whilst the effective 
action for the flavour modes scales as A^^ only if Np ~ Ns ~ 0{N), 

^5 = ^, ^F = ^. (2.9) 

It is worth repeating that in our notation, the theory has 2A^5 complex fundamental 
scalars, and Np Dirac fermions in the fundamental representation. When Np = Ns, 
the matter content in the flavour sector is that of Np hypermultiplets in the language 
of A/" = 2 supersymmetry. Similarly when ns = Q and np = A, the adjoint sector has 
the matter content of the A/" = 4 theory. 



3. Low temperature large-A^ transition 

At low temperatures TR <C 1, the adjoint modes are exponentially suppressed and 
we may ignore Sadj- The same is true for most of the flavour excitations except for 
those with energies comparable to the chemical potential ^. We take /i to be positive 
without loss of generahty and study the behaviour of the system as /x is increased 
from zero. 



- 8 - 



3.1 Truncated model: Light scalars ms <^ mp 

Taking tlie scalar modes to be strictly lighter than the fermions, for non-zero masses, 
we begin our low temperature analysis by keeping only the lightest mode carrying 
baryon number. The lightest scalar mode charged under U{1)b is the zero mode 
with £ = and energy Eq = y^SR~^~+^rn^, and we focus attention on the regime 



< /i ~£o = ySi?-2 + m2 , Ti?<l. (3.1) 
To be precise we work in the limit 

\R{fi - eo)\ ^ , TR^O with = fixed. (3.2) 

In this limit, the anti-squarks are Boltzmann suppressed by factors of e~^^^°~^^^ and 
can be ignored, as can all heavier fluctuations. It is then consistent to truncate the 
effective action to the lightest squark mode alone {U{1)b charge +1) so that 

Ses[U] ^ NusTr ln(l -(U), ( = e'^^^"-'^) , (3.3) 

with (, the effective fugacity. Within this truncation, the effective action describes all 
gauge-invariant states made from the lightest flavour modes {£ = scalar harmonics) 
on 5*^. In this regime it is clear that the effective matrix model for U has a complex 
action, requiring a non-standard approach towards its large- iV solution. 

3.2 Holomorphic matrix integral 

Rewriting the matrix integral for U in the eigenvalue basis ( |2.1|) , and including the 
Jacobian for this transformation - the well known Vandermonde determinant, we 
obtain, 

5cff[^a] = - J]ln sinf^^j +iVz.5 5^1n(l-Ce^'") + zQiV^e„. (3.4) 

ab=l ^ ^ a=l a=l 

Here we have explicitly implemented the condition 9a = 0, via a Lagrange mul- 
tiplier Q, enforcing the requirement that U has unit determinant. In the large-A^ 
limit, we expect that the integral over 9a is dominated by a saddle point. However, 
since the integrand is complex, we must allow for complex saddle points. In particu- 
lar, assuming that there exists such a saddle point which dominates the integral at 
large- iV, the corresponding equation of motion is, 

1 r0h-9a\ e''" 



a=l i^b) 



w ^ J+^^'C'r^^'G. (3.5) 
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The equation has solutions only if da are complex. To this end, we analytically 
continue the equation by introducing holomorphic variables, 

N 

Za = e'^- , with ^\^Za = . (3.6) 

0=1 

In terms of these variables the saddle point configuration satisfies 

1 ^ I 

J- V"^ Zb-\- Za Zb 

a=l {^b) 

Summing over 6, from this equation we learn that the Lagrange multiplier Q (times 
N) has the interpretation of average baryon number in the grand canonical ensemble, 

^ 1 C Za T dlnZ , ^ 

The squark number is N'^Q. For low temperatures, when // < £0, the fugacity is 
exponentially small and the squark number should be suppressed. As fi is increased, 
we expect something interesting to happen since there is the possibility of a pole 
in the above expression for Q leading to a divergent baryon number, signalling the 
onset of Bose condensation. 

3.3 Low fjL confined phase 

To solve the saddle point equation at large N, we take the eigenvalues Za to lie on a 
contour C in the complex z-plane. When the fugacity vanishes we know that C is the 
unit circle - this follows from the repulsive pairwise Vandermonde potential between 
the angular coordinates 9a- In the absence of any other potential terms at zero 
temperature and zero chemical potential, the eigenvalues spread uniformly around 
the unit circle. This is the confined phase of the theory. When the temperature is 
small, with a non-zero chemical potential and ^ < 1, we look for a closed contour C 
which is an appropriate deformation of the unit circle. 

It is natural to introduce a parameter t, with — tt < t < tt and treat the positions 
of the continuum of eigenvalues as a function z{t). When the fugacity C vanishes and 
C is the unit circle we must have z{t) — e**. More generally, we require -z(— tt) = z{Tr), 
for C to be a closed contour. In this continuum limit all discrete sTims can be replaced 
by integrals over C weighted by an eigenvalue density function p{z): 

(3.) 



c 



=1 

The eigenvalue density satisfies two constraints - a normalization condition (fixing 
the number of eigenvalues to N) , and the unit determinant condition on U : 
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The second requirement actually fixes the orientation of the branch cut of In z, so 
that it extends along the negative real axis from the origin to z = — oo. 

When the eigenvalues are uniformly distributed on the unit circle (the case with 
vanishing (), p{z) = ^. In this parametrization, the simple pole at the origin (en- 
closed by C) ensures that the eigenvalue density is normalized correctly. As the 
fugacity is increased smoothly from zero, as long as C remains closed, we do not ex- 
pect any further singularities to appear inside C. In the continuum limit the saddle 
point equation (|3l7| ) is 



r<pl^p{z')'-±^ + usCy^ = Q, zee, (3.11) 

' 2m z — z' 1 — Q z 



c 



where the integral over the contour C is to be treated as a principal value since the 
integrand has a singularity when z' = z along the contour of integration. 

It is fairly straightforward to see that, for small (, p{z) can be uniquely fixed to 

be. 



C small. (3.12) 



Even without an exphcit knowledge of C (which can now be deduced) , the assumption 
of analyticity has allowed us to determine p. 

The solution is valid as long as the simple pole at z = ^ lies outside C. This 
is certainly true when ( vanishes and should continue to be the correct picture as ( 
is increased smoothly. However, we expect something interesting to happen as ( is 
increased and the pole at z = ^, approaches the contour C on which the eigenvalues 
actually he. Note also that, as long as this pole lies outside C, the baryon number Q 
(defined in (|3.8| )) vanishes identically by Cauchy's theorem. 

Since p{z) is fixed, we can now determine C using Eq. ( |3.9|) , 



e'* = 7 -— , ~7T<t<7r (3.13) 

(l-C^)-s' - ^ ^ 

which implicitly fixes z as a function of t. As expected, when ( vanishes, we recover 
the unit circle z = e**. 

In this phase the expectation values of all the "holomorphic observables" , Tr ?7" 
vanish identically by the residue theorem, whereas the operators Tr U''^ are all non- 
vanishing. 



(3.14) 



This is a consequence of the difference between states carrying baryon and anti- 
baryon number (made from £ = modes of the scalars). The baryon number chemical 
potential biases the system towards baryons as opposed to anti-baryons. 
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3.4 Phase transitions 



The form of the contour C, and the location of the poles and zeroes of p{z) relative 
to it, signal the possible onset of phase transitions as we show below. For general 
values of i>s and (, p{z) has one simple pole and one simple zero in the 2;-plane, 

Pole : Zoo = 7 , Zero : Zq = ——^ r . (3.15) 

Interestingly, both lie on the real axis. Now, there are two qualitatively distinct 
possibilities depending on whether Zq is positive or negative. 




Figure 2: The relative positions of the contour C and the pole Z^o and zero Zq with 
increasing for vs < 1- The pole repels the eigenvalue distribution, and eventually the zero 
hits the contour C. Further increase in causes the distribution to split into disconnected 
components C and C" . 



3.4.1 2Ns<N: splitting eigenvalue "bubbles" 

Figure(0) shows the contour C, and the positions of the pole and zero as the fugacity 
(or the chemical potential) is increased. When z/5 < 1, Zq and Z^o lie to the right of 
C. However, as ( is increased both the pole and the zero approach C; the distribution 
of eigenvalues appears to be "repelled" by the pole, splitting the distribution in two 
disconnected contours when eventually the zero actually crosses the contour C. This 
splitting process leads to a non-analyticity in the eigenvalue density, suggestive of a 
phase transition. 

The critical value of ( for the onset of this behaviour can be obtained by substi- 
tuting Zq into the implicit expression ( p.l3[ ) for the contour. The requirement that 
C be real the leads to the solution, 

(1 _ ugys-i 

Ccrit = ^-T- , t = ±Us7l, 1^S<1- (3.16) 

Therefore there are two values of the parameter t (= ±2/57?), which simultaneously 
map to the point Zq when ( = Ccrit- Hence, two points on the contour C come 
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together at z = Zq, when the fugacity is tuned to the critical value. This explains 
the behaviour of the contour C encountered in Figure (0). In particular, 



C" = {z{t) ; —usn <t< usTi} 

C = {z{t) ; — TT < t < —UsTi , I's'^ < t < tt} 



(3.17) 



We infer from the extents of the two distributions as a function of t, that a fraction 
i^s of the eigenvalues {2Ns out of a total number N) have peeled off from the original 
contour C. 

For large (, Eg. (|3.13| ) has two qualitatively distinct solutions for z{t), 

z{t) =^ - 1 - — — , -T^vs <t<nus, (3.18) 



C^-'^s exp 



f-sgn(t)7ri/g 
I-US 



USTC < \t\ < 71 . 



The first of these has z ~ <^ 1, whilst the second has \z\ ~ (^'^s/i-'^s ^ x. 
These are the two contours C" and C in the limit of asymptotically large fugacities. 
The 2Ns eigenvalues making up the smaller contour C" have essentially collapsed 
to a point with parametrically small VEVs, whilst those lying on C are effectively 
uniformly distributed on a large circle. Correspondingly, the magnitude of the eigen- 
value density \p{z)\ exhibits a discontinuity for ( > Ccrit as shown in Figure (§). 
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Figure 3: The absolute value of p{z{t)) as a function of t as the fugacity is increased, 
for us = 0.5. When ^ = 2, the density vanishes at t = it^, whilst for larger values, a 
discontinuity develops, between the (high) density along the smaller contour C" and the 
larger one. The magnitude of the density along C remains small but always non- vanishing. 



The physical interpretation of this non-smooth behaviour in the eigenvalue dis- 
tribution is not immediately obvious. One possibility, which we have not explored, 
is that beyond ( = Ccrit a new saddle point solution with a gapped eigenvalue dis- 
tribution appears, on disjoint open contours. For the moment we will simply follow 
our original saddle point configuration beyond Ccrit, where it sphts, and attempt to 
understand its physical interpretation. 
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Despite the singular behaviour of the eigenvalue distribution, gauge invariant 
quantities such as the baryon number defined via Eq. (|3.8| ), (Tr?7") and {TrW"-) 
remain unchanged as ( is dialled past Ccrit- Here we are specifically referring to 
multiply wound loops with n fixed as is taken to infinity. The contributions to 
these quantities from the pole at z = cancel between the two disjoint contours. 
This is because C and C" must naturally be assigned opposite orientations when 
computing physical observables using the density function p{z). In particular, the 
total baryon/squark number remains zero during this process. The physical reason 
for this can be traced to the fact that the only matter modes we have kept in our 
discussion are the squark zero modes. When 2Ns < N, no gauge-invariant state 
carrying baryon number can be made solely using the complex scalars and therefore 
baryon number cannot change in this regime. We will see later how this picture is 
altered by the inclusion of additional degrees of freedom in the effective description. 

While the total baryon/squark number is zero, the two components of the eigen- 
value distribution each carry equal but opposite baryon number: 

NQ[C"] = -NQ[C'] = Nusil-us). (3.19) 

It is interesting to ask whether the splitting of the eigenvalue distribution into 
smaller "bubbles" , can be characterized by means of gauge- invariant order parame- 
ters as a phase transition. Possibly, Polyakov loops in rank k tensor representions, 
symmetric or antisymmetric, with k ~ 0{N), could exhibit specific non-analyticities 
associated to the splitting of the distribution. Such observables are known to encode 
gapping transitions of the Gross- Witten-Wadia (GWW) type [0, and can be com- 
puted using methods explicitly discussed in [^. We postpone further study of these 
issues for future work. 

Noting that Ccrit is strictly greater than unity for z/5 < 1, a surprising feature of 
the saddle point configuration above, is that it does not exhibit any obvious non- 
analyticities when C = 1- When C > 1; the chemical potential /i exceeds the energy eo 
of the lightest mode and, in a free theory, should destabilize the theory by inducing 
runaway VEVs for the squark zero modes. In particular the value of /i ~ £0 at which 
the splitting of the distributions occurs (at low temperatures), 

/i = £o + Tin Cent, (3.20) 

is strictly greater than Eq, and yet we did not encounter any sign of an instability as 
the /i was dialled past eo- 

One possible reason why the putative instability is not visible in the finite volume 
theory, has to do with gauge invariance. The squark fields are not gauge-invariant 
and when 2Ns < N, it is not possible to form gauge-singlet baryon operators in 
the matrix model truncated to the zero modes. Therefore no gauge-invariant order 
parameter exists to signal the onset of an instability (e.g. a runaway VEV) in the 
truncated model. 
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From the point of view of the effective matrix model for [/, a squark VEV 
instabihty is evaded by non-zero VEVs for the Polyakov loop eigenvalues in the 
confined phase. In the free theory, the action for the scalar zero modes, in the 
presence of the chemical potential and a diagonal VEV for Aq as in has the 

form 

Ns 

Cs = ^idrq] + iTa-qj + nqj) {drQj - iTa ■ qj - nqj) + {ml + ER-'^)qjqj + {q -> q) . 

(3.21) 

While the chemical potential acts as a negative mass squared for the scalar zero 
modes, non-vanishing eigenvalues of the Polyakov loop (a = {9a}) appear to want 
to stabilize the theory. 

3.4.2 Gross- Witten-Wadia transition for 2/5 > 1 

When the number of charged scalar flavours 2Ns is larger than A^, the number colors, 
we encounter a very different picture, with a new continuous deconfinement phase 
transition. Now the zero Zq, of the analytic function p{z), lies on the negative real 
axis. As ( is increased, Zq moves toward, and eventually sits on the contour C, 
precisely when the fugacity approaches the critical value 

_ {us - ir^-' 



c 



crit 



< 1 for z/5 > 1 . 



(3.22) 



The critical value follows from Eq. ( p.l3| ), upon substituting in the expression for 
Zq. Figure (^ illustrates this phenomenon for vs = 2. At this point the eigenvalue 
distribution has a zero and wants to develop a gap. Finding the new distribution 
and establishing its existence requires a separate analysis since the condition that C 
be closed, needs to be relaxed. Before we calculate the gapped distribution for the 
new phase, we explain the physical origin of this transition. 
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Figure 4: For 1/5 > 1 (number of squark flavours larger than the number of colors), the 
eigenvalue distribution develops a gap as (" is increased past Ccrit- The rightmost figure is 
just a continuation of the ungapped result beyond Ccrit- The correct gapped distribution is 
determined separately below. 



3.4.3 Baryons "melting" at large- AT^ 

The putative non-analyticities we have seen so far, in the large- A^ free theory in finite 
volume, can be provided a physical interpretation by looking at the grand-canonical 
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partition function of gauge-singlet states carrying baryon number. We first look at 
the case where the number of squark fiavours is larger than the number of colours'^. 
In this case (1/5 > 1), the theory has gauge-invariant excitations carrying baryon 
number, 

-D5(iii2...iiv) - (In 1i2 ■■■% ■■■liN e«ia2...aiv • 

The subscript '5" denotes that the baryon is made of scalars alone; in fact, from the 
zero modes of squark fields carrying baryon number In our notation these are 
the Qi and q\ with i = 1, 2, . . . Ng- In the truncated low temperature model, when 
2Ns > N, the degeneracy of single baryon states is 

^{Bs) = (^^'^ ^ew[N{uslnus-ius-l)ln{us-l))]. (3.24) 

Therefore the degeneracy of baryon-like states grows exponentially with N in the 
Veneziano limit. Each of these states has mass Nsq and U{1)b charge equal to N. 
Hence, the partition sum for gauge-invariant states at low temperatures wherein the 
"lightest states" charged under U{1)b contribute goes as, 

Z = 1 + exp [iV (z/5lnz/5 - (ug - 1) \n{iys - 1))] x 3-^'^^'°-^'^ + ... . (3.25) 

In the large N limit, the partition function diverges exponentially with N, when 

l^>eo-T{us In i^s - {i^s - 1) H^s - 1)) (3.26) 

which coincides with the critical value of the fugacity ( {us > 1) for the onset of a 
putative GWW transition in our matrix model. We will see below that this phase 
transition is similar to a deconfinement transition for squarks. At large N, due to 
the exponential growth of baryon operators, beyond the critical value of the chemical 
potential baryons are not the correct degrees of freedom, and the theory enters the 
phase of cold deconfined squarks. 

When the number of squark fiavours is smaller than the number of colours, it 
is not possible to form a baryon operator out of the squark fields (such operators 
vanish identically for 2Ns < N). However it is still possible to form adjoint-baryons, 
if we include adjoint modes. On the three-sphere, such exotic baryon-like states are 
possible even if the theory has no matter fields in the adjoint representation. The 
harmonics of the gauge field naturally provide such objects, and can be used to make 
adj oint-baryons . 

In particular, using a field $ transforming in the adjoint representation, we can 
make composite fields of the form and ^""qj, transforming in the fundamental 
representation of SU{N). It is easily seen that the lightest adjoint-baryon we can 
construct for 2A'^s < < ANs, is 

^The arguments presented here grew out of a stimulating discussion initiated by S. Minwalla. 
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For theories with 4A^5' < or 2z/s < 1 we will also need to include states made from 
contractions of higher powers of the adjoint mode, such as ^^^j, ^^Qi, etc. Taking 
the mass of the field $ to be e (this could be a scalar zero mode or the lightest 
mode of the gauge field), the lightest of the adjoint-baryon operators has a mass 
N{eo + (1 — z^s)^). The degeneracy of these adjoint baryons also grows exponentially 
with N as^. 



TV In 



(3.28) 



The exponential growth implies a phase transition when the chemical potential is 
increased beyond the critical value, 

/icrit = £o + (1 - i^s)e - ^ In S)(5Ad) . (3.29) 

It is interesting to see how the onset of this transition and the behaviour of the theory 
beyond it, is captured by the matrix model description. Of course, this will require 
going beyond the truncated model for the "lightest modes" in the presence of the 
chemical potential. In fact if the Kaluza-Klcin harmonics of the squarks are lighter 
than all adjoint modes, similar phase transitions can be triggered by the exponential 
growth of baryon states made from squarks carrying different angular momentum {£) 
quantum numbers. We postpone further discussion of this to Section 4. 

In order to better understand the nature of phase transitions induced by a chem- 
ical potential in this theory , we return to study the truncated matrix model in its 
gapped phase for 2Ns > N. 

3.4.4 Gapped Phase 

Let us now look for the gapped solution to the saddle point equation, which we expect 
to find when 2A^5 > and the fugacity exceeds the critical value beyond which the 
large- degeneracy of baryon operators leads to a divergent partition function. The 
contour C along which the eigenvalues are distributed, must be taken to be open. 
We begin by defining the resolvent function, 

U;(z) ^-iy P{Z') ^ . (3.30) 

The resolvent function will have a branch cut along the contour C, and importantly, 
branch points corresponding to the end-points of C. In particular, the saddle point 
equation, when expressed in terms of the resolvent, shows that the branch points 
must be of the square root type: 

-]-[u{z + e)+u{z-e)] = Q-vsC^^^ ^ C . (3.31) 

2 1 - 



^ At this stage, for simplicity, we are assuming that the number of such light adjoint fields is one. 
If the number of adjoint species is d, the degeneracy of adjoint baryons grows as {^^2Ns) 
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Here z ± e are two points infinitesimally close to the cut C, on either side of it. Al- 
though, generally speaking, branch cuts of complex functions can be moved around, 
in the present case there is one special choice of C which is determined uniquely by 
( p.9|) once the spectral density function p{z) known. The density function is also 
given by the discontinuity of the resolvent across its branch cut: 

^[uj{z + e) - uj{z - e)] = z p{z) , zeC. (3.32) 

This is a very useful relation, as it allows us to re-express the expectation values of 
matrix model observables as contour integrals enclosing the branch cut along C. In 
particular, line integrals along C involving the spectral density p{z) can be replaced 
by contour integrals involving the resolvent function: 

- — u{z) - . . . = f p{z) . . . , C enclosing C . (3.33) 

A-Ki z Jc 27rz 

From our analysis in the ungapped phase, we also expect the density function 
p{z) to be symmetric under reflection about the real axis, so that p{z) = p{z*). 
Therefore uj{z) must be of the general form 



ujiz) = f{z) ^{z-a){z~a*) + Z/5C - Q (3.34) 

where f{z) is an analytic function, which is completely determined by the require- 
ment that uj{z) be regular ai z = ^ and the asymptotics at infinity and the origin, 
namely, 

\im.uj{z) = 1, \im u(z) = —1. (3.35) 



We find that 



The final ingredient in the solution is the baryon number Q in the gapped phase. 
This is fixed by demanding the unit determinant condition on the matrix U, after 
expressing Eq.( p.lO| ) as an integral over a contour C enclosing C, 

^\nz^ = 0. (3.37) 

Since uj{z) itself has no singularities other than the branch cut enclosed by C, the 
latter contour can be deformed to enclose the other singularities of the integrand - 
the poles at z = and z = 00, and the branch cut of Inz which must be taken to 
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run along the negative real axis. Careful evaluation of the corresponding contour 
integrals yields^''' 



c 



:q + i 



iS+l 



(3.38) 



The transition to the gapped phase occuring at C = Ccrit is a third order phase 
transition, similar to the GWW model. The derivative of Q with respect to /i, which 
is the second derivative of the grand potential, is continuous across the gapping 
transition, 

dQ _ T2 dHnZ 



1 



lnQ + ln(Q + 1) +ln 



(3.39) 



At the transition point where ( = Ccrit and Q is vanishing, the second derivative 
of the grand potential also vanishes and is continuous, while the third derivative is 
discontinuous and diverges when ( > (crit (see Figure(|I)). 

We can now see in Figured^) what happens when the fugacity ( is increased 
beyond the critical value ^crit towards C = 1- As C approaches unity (equivalently 
the chemical potential approaches the energy Eq of the lightest mode), the baryon 
number Q diverges. This is the onset of Bose condensation of squarks. Furthermore 
as Q — )■ oo, the contour C shrinks to a point because the branch points ai z = a and 
z = a*, both approach the point 2; = 1 on the real axis. 
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Figure 5: Left : A depiction of the gapped distribution as the fugacity C is increased 
towards 1. It shrinks to the point z = 1 in the limit C ^ 1, when the pole Z^o collides 
with the eigenvalue distribution. Centre : The baryon number becomes non-zero when 
C > Ccrit and eventually diverges at ( = 1, signalling the onset of Bose condensation. 
Right: Continuity of second derivative of the grand potential (or first derivative of Q) 
across the third order transition. 



The third order transition we have seen above has a natural interpretation as 
deconfinement of squarks or fractionation of baryons. The "holomorphic" observables 

^The relative position of C with respect to the branch cut of Inz, plays an important role in 
obtaining the correct result. 

^This relation between ^ and Q suggests an underlying combinatorial origin, from Stirling's 
approximation, n\ ^ . 
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Tr[/" which vanished in the ungapped phase, are now non-zero. Equations ( 3.14 ) 
and ( p.33|) imply that expectation values of these are encoded in the Laurent series 
expansion of uj{z) at infinity, 

oo ^ 

u{z) = -1 - 2 5^(iTr [/")-. (3.40) 



n=l 



For example, 

which approaches unity as C — ^ 1 and Q ^ oo. Since the chemical potential breaks 
the symmetry under baryon charge conjugation, the charge conjugate set of observ- 
ables Tr f/''"", should have different expectation values. Once again. Equations ( p.l4| ) 
and ( |3.33| ) fix them in terms of the Taylor series expansion of w(-z) around the origin: 



UJ{Z} = . 

n=\ 



2 ^ (^Tr t/t") z-^ , (3.42) 



so that (;^Tr U'^) = distinct from its holomorphic counterpart. Since all the 

Polyakov loops (singly and multiply wound), have non- vanishing expectation values, 
the gapped distribution describes a deconfined phase of squarks ^ . 



It is interesting to note that in the limit C — ^ 1, at the onset of the Bose condensed 
phase, the branch cut C shrinks to the point z = 1, and the resolvent simplifies, 

\imu(z) = (3.43) 
C-^i 1 - z 

In this limit all single trace expectation values become equal, (^Trf/") = (^Trf/t") = 
1. This means that the angular variables {9a} all vanish and the gauge group is "un- 
broken". When the {6a} are distinct and non-zero, the gauge symmetry is broken to 
?7(l)^-i. 



4. Incorporating heavy modes on 

In the discussion above we saw the appearance of a third order deconfinement tran- 
sition in the general situation where squarks (scalars in the fundamental representa- 
tion) are the lightest modes. The non-analyticity of the grand potential as a function 
of baryonic chemical potential was demonstrated when the number of squark flavours 
{2Ns) was larger than the number of colors (iV). In this situation the non-analytic 
behaviour of the grand canonical partition function is attributed to the exponentially 
large degeneracy (~ e^) of baryon operators. 

When the number of squark flavours is small, i.e. 2Ns < N, the situation is 
somewhat subtle since baryon operators made from squark zero modes alone, do 
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not exist. In this case, we have aheady seen that the partition function truncated to 
squark zero modes (the "hghtest" degrees of freedom with U{1)b charge +1) displays 
no non-analyticities as the chemical potential is increased (despite "bubble splitting" 
of the eigenvalue distribution). However, baryon-like states, namely adjoint-baryons, 
do exist when then theory has additional adjoint modes which are neutral under 
U{1)b- In the Veneziano limit, the number of adjoint-baryons of a given mass also 
grows exponentially with N, possibly leading to a phase transition (|3.29|). Below we 



will demonstrate, using the approach of holomorphic matrix integrals, the existence 
of deconfinement transitions also when 2A^s' < N, when additional modes on 5*^ 
are included in the description. For sufficiently light squarks, another effect can 
become important: it is possible to form baryons using squarks with differing £ 
quantum numbers. These also exhibit an exponentially large degeneracy, and should 
be responsible for transitions in the large- iV matrix model. 

4.1 The effect of adjoint modes 

We now explore what happens when light adjoint bosonic fields neutral under U{1)b, 
are included in the effective theory for U. At the outset, we emphasize that such 
adjoint modes will appear in all gauge theories - they can be viewed either as arising 
from harmonics of free vector fields on or from additional adjoint matter fields (as 
in A/" = 2 supersymmetric theories). Including their lowest harmonics in our effective 
description at sufficiently low temperatures, the effective action for the unitary matrix 
U is 

Ses = NusTr ln(l -(U) -dJ^^T^^ ^"Tr t/^" (4.1) 



71=1 



where d is the number of such adjoint species and e is the associated energy. Here, 
we are implicitly assuming that there is always a range of parameters for which the 
anti-squarks, and other harmonics of fields on S*^ can be consistently ignored. We will 
eventually determine values of the chemical potential for which the modes retained 
above are the most important ones. 

The fugacity is defined as before, always with reference to the energy Eq of the 
lowest squark mode charged under IJ{\)b- For the lowest harmonics of vector fields 
on S*^, (i = 6 and e = 2/R, while for conformally coupled massless scalars, e = 1/R 
(see Table(0)). In what follows, for simplicity we set d = 1 (one adjoint species). 
At this juncture we impose no special restrictions on the chemical potential or the 
fugacity. Notice that in the low temperature limit, x is exponentially small and the 
adjoint contributions appear to be exponentially suppressed. Despite this, we will 
see that when the fugacity ( is cranked up to sufficiently large values, the adjoint 
modes drive a GWW deconfinement phase transition. 
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The large- equation of motion for the eigenvalues of U, following from the 
modified effective action ( [4 .11) , can be written in holomorphic form as before, 



dz' 



^z + z' ( z z' 
V + x 



z — z' \z' — X z z — X z' 



+ f^ = S. (4.2) 

1 - Cz 



z eC. 



At low temperatures a; <^ 1 and ^ ^ 1, the contour C is closed and all the poles 
in p{z) can be safely taken to lie outside C. Interestingly, this equation can also be 
solved exactly to yield p{z): 

(4.3) 

The infinite sum over n can be rewritten in terms of the q-Pochammer symbol as 
A 

dz 




(C^; a;)oo- The contour C itself is then given by z{t), where, using Eq.( p.9|) 



To understand the physics of this system systematically, working consistently at 
low temperatures (x <^ 1), let us begin by keeping only the first two terms (n = 0, 1) 
of the infinite sum, 

pW-^ + C..(^ + ^4^ + ...). (4.5) 

In the language of counting baryon-like states, this means we only count adjoint- 
baryons made from 2Ns squarks and (A^ — 2A^5') "adjoint-squarks" of the form $ ■ g^, 
each with a single adjoint matrix. 

Clearly, this goes one step beyond the truncated model of Section 3. As before, 
we can first locate the poles and zeros of this function. For small enough C, and x, the 
two simple poles at Z'^ = and Z^-* = (xC)~^ are far from the (closed) contour 
C which encloses the origin. Our approximate form of p{z) also contains two zeroes, 
which, for x <^ 1, are at 

C(l - Us) X C(l - 2Z/5) 

The first zero lies to the right of C when 2A^5 < A^ and we have already observed its 
effect in the context of the splitting of the eigenvalue distribution. We denote the 
critical value of C at which the splitting of eigenvalues occurs as Ceil- ^^i^ value, 
the zero Zq ^ hits the contour C. 

The second zero is new and has appeared solely due to the inclusion of the 
adjoint modes in our description. It is on the negative axis and to the left of C if 
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2Ns < N < ANs; else it lies on the positive real axis. Let us focus attention on the 

(2) 

case with 2Ns < N < ANg. As ( is increased, moves towards the contour C and 
eventually collides with it when (using Eg. (^.4])) 



c = c, 



(2) 
crit 



X 



(4.7) 



The corresponding value of the critical chemical potential, 



/^crit = (4.8) 

£o + (1 - i^s)e - T {us Hus) - {2vs - 1) H^^^s - 1) - (1 - i^s) ln(l - i^s)) , 

coincides precisely with the value deduced from the exponential growth of the de- 
generacy of adjoint-baryon states (Eg. ( p. 29]) ) for the case 2A^5 < < 4^"^. The 
phenomenon is summarized in Figure (0) 



/ V = 2/3' 
1 (= 1-50 , 






mt 9 



V = 21 
f = 2.8f 


3^ 


9-X 









2/3 



f = 2.90 



Figure 6: Plots demonstrating the behaviour of the contour C of eigenvalues in the z-plane, 
with increasing fugacity and the relative postions of poles (in blue) and zeroes (in red) of 
p{z). At a critical value of after the "splitting" of the contour, the distribution develops 
a zero, signalling the onset of a gap. The curves were obtained assuming x = 0.1. 

The preceding analysis demonstrates that even when the number of sguark 
flavours is smaller than the number of colors, the theory on is subject to a low tem- 
perature GWW deconfinement transition when the chemical potential approaches a 
threshold value determined by the mass of the lightest adjoint-baryons and their de- 
generacy in the Veneziano limit. We also deduce that the effective action retaining 
only sguark zero modes and the adjoints is a good description at low temperatures 
when iV/i is in the vicinity of the mass of the lightest adjoint-baryon. 

4.2 Effect of higher squark harmonics 

When 2Ns < it is obvious that baryons cannot be made from sguark zero modes. 
However one may still use sguark modes with differing angular momentum guantum 
numbers to form gauge-singlet baryons. Such baryons could be formed by 2A^5 sguark 
zero modes {^ = 0) and A" — 2A^5 sguark modes with i = 1. Given a degeneracy d 
for the i = 1 scalar spherical harmonics, the number of such baryons is {]^^2Ns) ^^'^ 
much in the same way as for adjoint baryons, it should lead to a phase transition 
with increasing fugacity. 
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The baryons made from higher harmonics of squark modes become relevant only 
if the energies of the squark modes are significantly smaller than R^^ so that all 
adjoint harmonics (including the gauge fields) are heavier. In this situation the 
effective action for U, keeping only the £ = and £ = 1 squark modes, is 

S^s ^ NusTr \n(l -CU) + dN Us Tr \n(l-xCU) (4.9) 

^ — g-/3(£0-/i) J. — g-/3(£-£o) 

Here e is the energy of the i = 1 harmonic, and the degeneracy of this mode = 4. 
Now the saddle point configuration in the ungapped, confined phase, is solved by the 
density 

p{z) = - + — + — . 4.10) 

z i- — Qz 1 — xQz 

The density function has zeroes at 

In the range INs < N < {d+ l)2Ns, the first zero lies on the positive real axis, whilst 
the second zero is on the negative real axis, to the left of the eigenvalue distribution. 
With increasing chemical potential, a zero in the eigenvalue distribution first appears 
when, 

or 

= -.^0 + (1 - . - T In ( (..^^^..^.^/.'^^^^^ _,)„,.,„,_. ) . (4.13) 

The critical value of the chemical potential is determined by the mass and degeneracy 
of the lightest baryon operators composed of the squark modes with £ = and £ = 1. 
Below we will work out what happens to the system as the chemical potential is 
dialled past this critical value. In particular we will show the emergence of a gapped 
"deconfined" phase. 

4.3 Gapped Phase 

The transition to the gapped phase for z/^ < 1 (small number of flavours) is always 
preceded by an eigenvalue splitting transition. The eigenvalues which have split off 
and formed the smaller "bubbles" do not appear to play any role in the subsequent 
dynamics as the chemical potential is increased to large values. 
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4.3.1 Gapping/Bose condensation via higher harmonics 

Let us attempt to find the gapped distribution which results when the chemical 
potential is dialled past the critical value in Eq.( [4.12| ). Since the distribution has 



already split prior to the gapping transition, we write the spectral density function 
as a sum of two (disjoint) pieces, 

p{z) = pi{z), zeC (4.14) 

= P2{z) , Z eC" . 

where the first term is associated to the larger contour C, and the second term yields 
the density of the smaller "bubble" C" (see e.g. Fig.(|])). In the splitting transition, 
a fraction z/5 of the N eigenvalues break off to form C" . We must therefore have 

dz f dz 

From our earlier analysis, we know that the eigenvalue bubble is parametrically small, 
centred around z ~ foi^ C ^ 1- The configuration at the large- iV saddle point 
solves 

„ [ dz' , ,^ z + z' VsC z uqxCz ^ 

r — p^' 7 + + ^ 1 V = Q- 4.16 

JciijQi, ZTTt Z — Z 1 — C-2 I — XQZ 

We further recall from the discussion on the eigenvalue splitting transition that the 
contour C" must be traversed in the clockwise sense. In the low temperature limit, 
X <^ 1, with C, large ^, the spectral densities in the two disjoint contours can be 
determined as follows. For z G C, due to the small radius of C" (~ C~^~'^^)) the 
saddle point equation reduces to, 

^'/|^^P,M^ + <if^ = S (4.17) 

Jqi Zm z — z 1 — z 

where we have used the normalization condition that P2{z)-^^ = vs- Remarkably, 
this is exactly the equation we solved to find the gapped phase of the truncated 
model in Section p.4.4| , provided we make the following replacements: 



J^s^T^ Q^T^ C^xC. (4.18) 

1 - z/5 I-VS 

The factor of (1 — z/5) in the denominator accounts for the fact that pi{z) is a density 
for (1 — vs)N eigenvalues making up C . The resolvent function for the gapped 



^Note that C, — cxp(— /3(eo — /i)) where £0 is the energy of the ^ = mode. To see the gapping 
phase transition /i needs to be close to the mass of the baryon made from i — and 1=1 harmonics. 
Due to the finite energy gap between the ^ = and 1 = 1 modes, C, will naturally be parametrically 
large for low enough temperatures. 
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distribution pi{z) is 



1 r dz' , z + z' 
^ ~ t^s Jc' 27ri z — z' 



co.iz) ^ / — Pi(/)r^, (4.19) 



u,{z) = A(z)v/(.-al)(.-aO-^+ ^_^^^_^^^ . 

The function fi{z) and the positions of the branch points at ai and are given by 
Eq. (|3.36|) with the replacements indicated in Eq.( |4.18D . It is useful to verify that 



right after the onset of the gapping transition where Q — ?• 0, the two branch-points 
coalesce and 

lim ai = ^ , (4.20) 

coinciding with the position of the zero (see Eg. ( |4.11[ )) of the ungapped distribution 
at the onset of the phase transition. Therefore the solution we have found is indeed 
connected, across a continuous transition, to the ungapped/confined phase. 

To find the distribution p2, we pick a point z G C" and noting that on this 
contour z ~ the saddle point equation becomes 

f dz' , ,^ z + z' , ^ vqCz , ^ 

f — P2{z') J - (1 - us) + = . (4.21) 

Jc" 2711 z — z' 1 — z 

Here we have used Jc Pi{^)'i^ = (1 ~ '^s)- This equation is actually solved by the 
same distribution function we found in the ungapped phase of the model truncated 
to zero modes, 

A final, important consistency check of our picture follows from the relation 
between the baryon number (Q) and the fugacity (, in this phase. As in the ba- 
sic truncated setup, the relation can be derived by enforcing the unit determinant 
condition on U, 

dz 

— p{z)\nz = 0, (4.23) 



c'uc" 



and we find 




(4.24) 



This means that the baryon number diverges Q — > oo when (x^~'^^ = 1 or p = 
^s^o + (1 ~ ^s) ^- Therefore when the baryon number chemical potential Np equals 
the mass of the lightest baryon state that can be formed from the spherical harmonics 
of the squark modes (£ = and £ = 1 in this case), Bose condensation of squarks 
occurs. We believe that it is the squarks which condense, because slightly before this 
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happens (at low non-zero temperatures), the system enters a deconfined phase by 
passing through a third order GWW transition. 

In all the GWW transitions we have encountered, the baryon number A^Q is zero 
below the transition and begins to increase smoothly above it. Therefore, using Q = 
W^^§ir-> grand potential is a constant independent of /i below the critical value 
of the chemical potential, and becomes a non-trivial function of /i above it. It follows 
from our results, such as, Eq.( 4.24j) that the grand potential above the transition is a 



complicated function of the form — /3 \nZ = N"^ F{vs-, /i, T). In contrast, in the low ^ 
confined phase, we will see below that the grand potential is of the form ~ Ng F{T, /i), 
indicating that the contributions are from a gas of gauge singlet mesonic states. 

4.3.2 Gapping due to adjoint modes 

We have already shown that when the number of squark flavours is smaller than 
the number of colours, it is possible to make gauge-singlet adjoint baryons. The 
exponentially large degeneracy of these states drives the system towards a GWW 
"gapping" transition (see Section |0|). Below we sketch the steps involved in finding 



the gapped phase of the model when light adjoint modes are included in the effective 
description. Since we expect a finite energy gap between the adjoint states and 
the squark zero modes, the GWW transition is preceded by the breaking up of the 
eigenvalue distribution into disjoint contours C and C". As before we write pi{z), 
the gapped spectral density along C, in terms of a resolvent. 



f dz 1 f dz' z + z' 

Jq, 2m 1 — z^s Jc' z — z' 



(4.25) 



The resolvent, by definition, is analytic on the z-plane with its only singularity being 
a branch cut along C. The saddle point condition for the effective model with adjoint 
degrees of freedom, can then be rewritten as 

1 1 1 1 Z/Q 

--{ui{z + e)+ui{z-e)) + -uji{xz) + -uji{x-^z) = f ^, (4.26) 

Z Z Z i — Ug 

zeC , X = e-^' . 

We have used the fact that ~ —us for large (, and for z G C . This is an 

interesting condition as it can be recast in a different form by defining a new analytic 
function, 

G(.)..,(v/;.-)-.,(^)-2f±^j;il. (4.27) 

G{z) has two branch cuts (ignoring for the moment, the logarithmic branch point) 
and the saddle point equation can be expressed as a gluing condition across these 
two branch cuts: 

G{V^z ± e) = T e) . (4.28) 
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In principle, this condition and the asymptotics of G{z) uniquely determine the 
resolvent. We will not pursue the analysis further in this paper. Curiously, such 
a gluing condition defines a Riemann surface of genus one; any significance of this 
to the physical problem at hand is unclear. Very similar "gluing problems" arise as 



solutions to large- iV holomorphic matrix models in very different contexts p2 |. 



4.4 Fermion flavour modes 

We now turn to the physics of quarks or fermionic flavour modes on S^. The matrix 
model incorporating these at finite chemical potential has already been studied in 
detail in [|^. Technically speaking, the saddle point equations for the fermionic 
flavours are related to the bosonic ones by the replacement ( — ?■ —( and us — > —up- 
However, the phase structure is very different, as must be expected for states with 
different statistics. 

In supersymmetric theories on 5*^, fermions are heavier than their bosonic coun- 
terparts (see Table(|T]) in the Appendix). However, if the quarks happen to be the 
lightest modes charged under U{1)b with energy cq and degeneracy d, the number 
of (lightest) baryon states would be 



3(5) 



2Nf d + N 
N 



exp 



N In 



'i+duFy-^'^"^ 



d Up 



(4.29) 



Here Np is the number of Dirac fermion flavours. Now gauge-singlet baryon operators 
exist for any non-zero value of Np with degeneracies scaling exponentially with N . 
Consequently, we expect a phase transition when n = Sq — ^\n. T){B). Such a phase 
transition was indeed seen in In fact, the fermionic system displays two such 
transitions, one on either side of /i = eo at 



^crit 



£oTi In S)(5). 



(4.30) 



Both these transitions are third order. The first is a deconfinement transition of 
the kind we have discussed above, where the baryons "dissolve" into quarks and the 
distribution of eigenvalues of the Polyakov loop matrix acquires a gap. The second 
transition is less intuitive, but can be understood as a complete occupation of the 
lowest Fermi level following which the theory re-enters a confined phase. Across this 
transition the eigenvalues of the Polyakov loop matrix return from the gapped to an 
ungapped continuum distribution. 

Let us briefly re-derive the above results for fermion flavours, by making use of 
the replacement C ~^ ~C and vs in Eq. (p3|) . Assuming that the degeneracy 

of the lightest quark modes (i = 1, we follow the analysis in Section 3.4. The crucial 
difference with respect to the bosonic case is that the pole (at Z^o = — 1/C) and zero 
(at Zq = — 1/C(1 + i^p)) are always on opposite sides of the eigenvalue distribution in 
the ungapped phase. In addition, the zero always hits the contour C with increasing 
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( irrespective of the value of up, thus triggering a GWW transition. This is as 
expected for fermions, since the corresponding baryon operators Bp exist for any 
I'F- Similarly, adapting the formula for baryon number Eg. ( |3.38D to fermions, we 
find that Q increases from Q = to a maximum of Q = up a.s ( is increased from 
(-.^ = z/^^/(l + i/^)i+'-F to = (l + upf+'^^/iyp''. Finally we turn to the formula for 
the end-points of the gapped distribution ( p. 361 ) and modify it to suit the fermionic 
case. We immediately obtain that the the contour is closed, when Q = and also 
when Q = Up. These are the two successive GWW transitions. 

The main physical effect is that following the two third-order transitions, given 
a degeneracy d for the lightest fermion harmonic, baryon number jumps by d vp N . 
Such jumps (each accompanied by a pair of third order transitions) occur with in- 
creasing /i as successive fermionic energy levels are encountered and filled. We refer 
the reader to [E4l for further details. 



5. Finite temperature, small ^ regime 

We have explored the low temperature regime of the phase diagram in the /i — T 
plane. The remaining question is what happens to the system as the temperature 
is increased for any fixed yU. We know that the low temperature phase boundary 
between confined (or mesonic) and deconfined phases is a straight line whose slope 
is determined by the logarithm of the large-A^ degeneracy of the (light) baryonic 
states. For small /x, the finite temperature behaviour can be understood along the 
lines of [jll], |19|, ^ |2T|. Defining the "moments" of the eigenvalue distribution (the 
normalized Polyakov loops) 

p.^^Tr^", p_„^lTrf/tn^ (5.1) 

the effective action is 

oo 

5efr[p] = N^Y.- [{^-zv{x^)-nsZs{x^) + {-irnpZp{x^))p^p.r, (5.2) 

n=l 

- {vsZs{nl^. msR) - (-1)" Z^(n|, mpR)) (p„ e"'^'^ + e""^'^ p_„)] 

The free gauge theory with no flavours {vs = vp = ^) exhibits a first order de- 
confinement transition with increasing temperature In pure Yang-Mills theory 
{ns = Up = 0), this happens when zv{x) = 1, and the "lightest mode" pi is forced to 
condense. From this condition, the corresponding Hagedorn/deconfinement tempera- 
ture is TR = 1/ ln(2 + y3). Similarly for AT = 4 SYM {ns = Q,np = 4), the first mo- 
ment pi condenses when zv{x) + 6 zs{x) +Azp{x) = 1 or when TR= l/ln(7 + 4y3). 

In the theory with fundamental flavours the situation is slightly different. For 
generic non-zero p, both p„ and p_„ are non-zero and unequal, and their values 



- 29 - 



determined by solving the (linear) equations of motion, 



^ g-n/3M Zsjnj^, msR) - up (-1)" ^F(n^, mpR) 

for n 7^ and po = 1- The first point to note is that upon plugging the solution back 
into the action, as long as this saddle point is valid, we expect that the free energy 
is the sum of terms of (9(A^|), 0{Np) and 0{NfNs). This is indeed the case in the 
low temperature ungapped phase |T^, leading to the interpretation that the grand 
potential should be attributed to colour-singlet mesons. This is the confined phase. 

The onset of the deconfinement transition is deduced by locating (for any given 
/i), the value of the temperature at which a zero appears on the contour C of eigen- 
values. It is a complicated excercise in general, but somewhat easier for vanishing p. 
Setting /i = and examining the effective potential for the "lightest modes" pi and 
p_i (in an approximation which ignores all higher moments), we obtain, 

_ _ us Zsij, msR) + Zpij, mpR) 

Pi - p-i - -j 7—^ ■ ^^-^^ 

1 — zv[x) — ns zs[x) — np zf{x) 

The non- uniform distribution develops a gap when pi = p_i = i, so that the spectral 
density on the circle then takes the form p{6) = ^{1 + cos^^) with a zero at = tt. 
This is a third order GWW transition ||19|, and is the natural continuation of the 



low temperature phase boundary to finite temperature and vanishing p. 



6. The phase diagram: discussion and conclusions 

The general picture is now fairly clear: at low enough temperature, with increas- 
ing baryon number chemical potential, the theory on is subject to a third order 
GWW deconfinement transition when the baryon number chemical potential (A^p) 
approaches the threshold set by the mass of the lightest baryon state and its expo- 
nentially large degeneracy in the Veneziano limit. Depending on whether the theory 
has fundamental fermions or not, there can be more than one such deconfinement 
transition. When the mass of the lightest scalar baryon (a baryonic operator made 
purely from fundamental and adjoint scalars) is approached, a third order baryon 
melting transition ensues followed by the onset of Bose condensation of the squark 
degrees of freedom. The free gauge theory is not defined beyond this point. We have 
also shown that in a different region of the p — T plane, namely, low p and finite T, 
the theory undergoes a third order transition as a function of temperature. Putting 
these pieces together we arrive at the general form of the phase diagram depicted in 
Figures (|1|) and (0). 

In both figures, the shaded region with Nfi > is a region of instability for 
the free theory, where the grand canonical ensemble is ill-defined (e.g. O, IT^)- It is 
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Figure 7: The phase diagram for the free theory, assuming that the hghtest baryon states 
are constituted of squarks. (See Figure(l) for the general situation when quarks are hght.) 



conceivable that upon introducing appropriate self-interaction potentials, the shaded 
region may have a stable ground state with a VEV for a gauge invariant operator, 
breaking U{1)b spontaneously. It can be shown that for SUSY QCD theories or 
the theory corresponding to the D3-D7 system (A^ = 2 SYM with hypermultiplets), 
there are flat directions in the scalar potential (Higgs branches and mixed Coulomb- 
Higgs branches) which are then responsible for runaway unstable directions when fi is 
non-zero. For any non-zero temperature and small coupling however, at the origin of 
the field space, the squarks and adjoint scalars will obtain a temperature dependent 
thermal mass so that ms{TY ~ fn^s + ^^{T)-, where A is the 't Hooft coupling and 
5{T) is exponentially small (~ e^^/™) at low temperature [|T3| and quadratic in T 
at high temperature. This suggests that, with increasing temperature, the chemical 
potential will need to be increased to a higher threshold value before a runaway 
instability kicks in. Therefore we may expect a small region in the fi — T plane, 
bounded by Nfi = tjib and the dashed line in Figures(|I]) and (^, where a metastable, 
deconfined phase exists and where all scalars have vanishing expectation values (see 



1^, |T5| for analogous discussions in the case of TV = 4 SYM with R-symmetry 
chemical potentials). 

At strong coupling (and in flat space), the unstable region has been studied in 
the context of the D3-D7 system for a small number of flavours and it should 
be possible to investigate this within the framework of which incorporates the 
back-reaction from flavour branes. In the 'probe' approximation (fixed number of 
flavours in the large- limit), the phase diagram in the n — T plane of the D3-D7 
system was extensively explored in Despite being in the probe hmit. 



we would expect a Bose condensed phase (in addition to a Fermi surface from the 
quarks) due to squark modes charged under U{1)b in the setup (see |^ ^ ^ for 
related discussions). An extensive search for instabilities and U{1)b breaking in the 
world- volume fluctuations of a probe D7-brane in AdS^ x has recently been carried 



out in M], with a negative result. It would be interesting to understand this better 
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within a general setup of the kind we have studied but with (small) interactions 
switched on. It would also be of interest to find what role, if any, is played by states 
such as adjoint-baryons in this picture. 

The low temperature, finite density transitions we have found in the Veneziano 
large- iV limit, are due to the exponentially large degeneracy of baryon states. It is 
reasonable to expect these to persist for any coupling, at large- iV (e.g 0). If so, 
within an appropriate dual gravity set-up, this would correspond to a black hole 
solution dual to the phase of cold, dense, deconfined (s) quarks. 

A final comment should be made about the phases that we have found (Section 
4.2) with order parameters that involve modes of the fields that carry non-trivial 
angular momentum . These inhomogeneous phases should involve the spontaneous 
breaking of rotational symmetry on the S"^. This feature would be worthy of future 
study. 
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Appendix A: Harmonics on 



Table 1: Harmonics on the three-sphere 



Field 


Angular momentum {£) 


Energy (ee) 


Degeneracy (di) 




e>o 




2£{i + 2) 


a 


i>0 




{i + iy 


Ao 


i>0 






Fermions 


e>o 




2i{i+l) 


Massless scalars 




i£+l)R-^ 




(conformal coupling) 









We list the energies of Kaluza-Klein harmonics of scalar, vector and fermion 
fields on a spatial three-sphere. For details of the spherical harmonic decomposition 
we refer the reader to the discussion in . The main point to note is that fermions 
ip, and the spatial modes of the vector field denoted by Ai {i = 1, 2, 3), do not have 
zero modes on S^. In general, the spatial components of a gauge field Ai can be 
split into the image and kernel of Vj (the former are the transverse, physical degrees 
of freedom). In particular A^ = + C* with C* = V*/ and Vj-B* = 0. The time 
component of the gauge field A^ and scalar fields are allowed to have zero modes on 
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